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Abstract 

A model-independent analysis of the infinite-momentum-frame charge density of partons in the 
transverse plane is presented for the nucleon. We find that the neutron parton charge density is 
negative at the center, so that the square of the transverse charge radius is positive, in contrast with 
many expectations. Additionally, the proton's central u quark charge density is larger than that of the 
d quark by about 70%. The proton (neutron) charge density has a long range positively (negatively) 
charged component. 
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A truly impressive level of experimental technique, effort and ingenuity has been brought 
to measuring the electromagnetic form factors of the proton and neutron (nucleon) [lj. These 
quantities are probability amplitudes that the nucleon can absorb a given amount of momentum 
and remain in the ground state, and as such should be important sources of information about 
the nucleon charge and magnetization densities. 

The text-book interpretation of these form factors is that their Fourier transforms are mea- 
surements of the charge and magnetization densities. But the initial and final nucleons have 
different momentum, and therefore different wave functions. This is because the relativistic 
boost operator that transforms a nucleon at rest into a moving one changes the wave function 
in a manner that depends on the momentum of the nucleon. The presence of different wave 
functions of the initial and final nucleons invalidates a probability or density interpretation. 

A proper determination of a charge density requires that the quantity be related to the 
square of a wave function or of a field operator. The technical solution to the problem of 
determining the relevant density operator has been known for a long time 0], and has been 
elegantly explained recently [1,0] in terms of generalized parton distributions GPDs. The charge 
density p(b) \^ of partons in the transverse plane is a two-dimensional Fourier transform of 
the Fi form factor. Here we present the first phenomenological analysis of existing data to 
determine p(h) for the neutron and proton. The results for the neutron contradict the long- 
standing notion, derived from both gluon-exchange and meson-cloud models 0, 0], that the 
non-vanishing charge density at the center of the neutron is positive. 

We begin by by presenting definitions of the form factors. Let J M (x) be the electromagnetic 
current operator, in units of the proton charge. Then the nucleon form factors are given by 

(p>,\'mO)\p,\)=u(p',\') ^F 1 (Q 2 )+ l —q a F 2 (Q 2 ))u(p,\), (1) 

where the momentum transfer q a = p' a — p a is taken as space-like, so that Q 2 = —q 2 > 0. 
The nucleon polarization states are chosen to be those of definite light-cone helicities A, A' [§]. 
The charge (Dirac) form factor is F±, normalized such that Fi(0) is the nucleon charge, and 
the magnetic (Pauli) form factor is F 2 , normalized such that F 2 (0) is the anomalous magnetic 
moment, the Sachs form factors [91 



G E (Q 2 ) = F^Q 2 ) - J^ 2 (Q 2 ), G M (Q 2 ) = F^Q 2 ) + F 2 (Q 2 ), (2) 

were introduced so as to provide an expression for the electron-nucleon cross section (in the 
one photon exchange approximation) that depends on the quantities G\ and G 2 M but not 
the product Ge Gm- In the Breit frame, in which p = — p', Ge is the nucleon helicity flip 
matrix element of J°. Furthermore, the scattering of neutrons from the electron cloud of atoms 
measures the derivative —dGE(Q 2 )/dQ 2 at Q 2 = 0, widely interpreted as six times the mean- 
square charge radius of the neutron. However, any probability or density interpretation of Ge 
is spoiled by a non-zero value of Q 2 , no matter how small [lOj. Any attempt to analytically 
incorporate relativistic corrections in a p 2 /m 2 type of expansion would be doomed, by the 
presence of the quark mass, m q , to be model-dependent. 

The newly appreciated and widely studied generalized parton distributions GPDs are 
of high current interest because they can be related to the total angular momentum carried 
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by quarks in the nucleon and can be determined using deeply virtual Compton scattering 
experiments IjJ. These distributions are specific matrix elements of quark field operators, 



between nucleon states, which in contrast to the usual quark distribution functions, do not 
have the same momenta. We consider the specific case in which the longitudinal momentum 
transfer £ is zero, and A' = A. Then, in th e lig ht- cone gauge, A + = 0, the matrix element 
defining the GPD, H q for a quark of flavor q |l3|] is 

H q (x,t) = /^(p + ,p',A|g(-^,0) 7 + g(^,0)b + ,p,A)e^. (3) 

We use the abbreviation H q {x,^ = 0,t) = H q (x,t) and —t = —{p' — p) 2 — (p' — p) 2 = 
—q 2 = Q 2 . The simple form of t results from its invariance under transverse boosts [8|: Lorentz 
transformations, defined by a transverse vector v that transform a four-vector k according to 
k + — > k + , k — > k — k + w and k~ such that k 2 is unchanged. These quantities are part of a 
kinematic subgroup of the Poincare group that obey the same commutation relations as those 
among the generators of the Galilean transformations for non-relativistic quantum mechanics in 
the transverse plane. The presence of the operator 7+ insures that independent field operators 
appear in the matrix element. 

GPDs allow for a unified description of a number of hadronic properties [l2{]. The most 
relevant for us are that if t = they reduce to conventional PDFs H g (x,0) = q(x), and that 
the integration of H q over x yields the nucleon electromagnetic form factor: 

= E e <? / dxH q {x,t). (4) 



The spatial structure of a nucleon can be examined if one uses 0, 0, H[ nucleonic states that 
are transversely localized. The state with transverse center of mass R set to is formed by 
taking a linear superposition of states of transverse momentum. In particular, 



R = o,a) = at/^|p + ,p,a). (5) 



where \p + , p, A) are light-cone helicity eigenstates [8( and M is a normalization factor satisfying 
|A/f / Q$ = 1. References 0, M use wave packet treatments that avoid states normalized 



(27T) 2 

to S functions, but this leads to the same results as using Eq. ([HJ). Note however, the relevant 
range of integration in Eq. (jSJ) must be restricted to |p| p + to maintain the interpretation 
of a nucleon moving with well-defined longitudinal momentum [l~4| . Thus we use a frame with 
very large p + . It is in just such a frame that the interpretation of a nucleon as a set of a large 
number of partons is valid. 

Using Eq. (jSJ) sets the transverse center of momentum of a state of total very large momentum 
p + to zero, so that transverse distance b relative to R. can be defined. Thus we may generalize 
the quark-field operator appearing in Eq. ([3]) by making a translation: 



6 q (x,b)= /^ g t_ * b ?+ V. b ) e ■ ^ 




3 



The impact parameter dependent PDF is defined [l4j as the matrix element of this operator in 
the state of Eq. (J3J): 



q(x, b) = (p + , R = 0, A 6 q (x, b) p + , R = 0, A) . 



(7) 



The use of Eq. (jSJ) in Eq. (0) allows one to show that q(x, b) is the two-dimensional Fourier 

d 2 q 



transform of the GPD H q : 



?(*,b) = / -^e^ h H q (x,t = - q 2 



(27T) 

with i/g appearing because the initial and final helicities are each A. A complete determination 
of H q (x,t) (with t < 0) would determine g(jr, b). 

One finds a probability interpretation [2[] by integrating b) over all values of x. This 
sets the value of x~ to 0, so that 



J dx q(x, b) = R = 0, A| q\_(Q, b)g+(0, b) R = 0, A) . 



(9) 



If one multiplies the above relation by the quark charge e q (in units of e), sums over quark 
flavors, uses Eq. (j^D with O q (x, b) = e~ l ^' h O q (x, 0)e l6b along with Eq. (jlj), the resulting infinite- 
momentum-frame IMF parton charge density in transverse space is 



d 2 



(2tt)* 



q 2 )e i qb . 



(10) 



This transverse charge density p(b), based on iq {2J, seems to contrast with the lore relating 
the charge density to Ge- There is no conflict, a feature understood by considering the relation 
between the Breit frame and the IMF. In the Breit frame, the helicity flip matrix element of 
J° is Ge- One uses a Lorentz transformation to consider the form factor in the IMF. Then the 
operator helicity-flip matrix element of J° becomes the helicity non-flip matrix element of J + , 
or Fx- The transformation to the IMF gains model independence while replacing Ge by F±. 

We exploit E q. (fTOl) by using measured form factors to determine p(b). Recent parameter- 



izations [16|, [17|, [18| of Ge and Gm are very useful so we use Eq. (j5J) to obtain F\ in terms of 



Ge, Gm- Then p(b) can be expressed as a simple integral of known functions: 



P(b) 



dQQ - UQ b)^ 



tG m (Q* 



2tt 



1 + r 



'11^ 



with r = ^2 and Jo a cylindrical Bessel function. 

A straightforward application of Eq. ( TlTT) to the proton using the parameterizations 16|, ll7 
yields the results shown in the upper panel of Fig. 1. The curves obtained using the two different 
parameterizations overlap. Furthermore, there is negligible sensitivity to form factors at very 
high values of Q 2 that are currently unmeasured. The density is peaked at low values of b, but 
has a long positive tail, suggestive of a long-ranged, positively charged pion cloud. 

The neutron results are shown in the lower panel of Fig. 1. The curves obtained using the 
two different parameterizations seem to overlap, but see below. The surprising result is that 
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FIG. 1: (Color online) Upper panel: proton charge d^si^Fp|&).Lower panel: neutron charge density. 
The solid curves use the parameterization of 171 ]. and the dashed (red) curve uses [16]. 



the central neutron charge density is negative. If the neutron is sometimes a proton surrounded 
by a negatively charged pionic cloud, one would expect to obtain a positive central density 
Another mechanism involving correlations in the nucleonic wave function induced by one gluon 
exchange would also lead to a positive central density because the interaction between two 
identical d quarks^ is repulsive. The values of the integral of Eq. (fTTl) are somewhat sensitive 
to the regime 2 < r < 4 for which Ge is as yet unmeasured. About 30% of the value of p(0) 
arises from this region. 

The negative central density deserves further explanation. The upper panel of Fig. 2 shows 
Fi for the neutron obtained using the two different parameterizations which are observably 
different. However, in both cases F\ is negative (because of the dominance of the Gm term 
of Eq. (HIJ) for all values of Q 2 . This along with taking 6 = 0, J (Qb) = 1 in Eq. (TTT]) leads 
immediately to the central negative result. The long range structure of the charge density is 
captured by displaying the quantity bp(b) in the lower panel of Fig. 2. At very large distances 
from the center, again suggesting the existence of the long-ranged pion cloud. 

The present analysis provides detailed information about the location of charge density 
within the nucleon, and also incorporates the lore regarding mean-square-radii, MSR. It has 
long been known that the MSR defined by the form factor Ge is dominated by the Foldy term 
-1.91/(4M 2 ) = -0.126 fm 2 [j!, M ■ arising from the neutron magnetic moment ^(0). The 



experimental value of the Ge MSR, cited in [17J, is (-0.114 ± 0.003) fm 2 , so the MSR associated 
with F\ (obtained from the integral / d 2 bb 2 p{b)) is small and positive (+0.012 fm 2 ). This result 
is consistent with Figs. 1 and 2. However, knowing the MRS of F\ does not, by itself, allow 
one to conclude that the central neutron charge density is negative, does not reveal the critical 
model-independent feature that at the very largest distances the charge density is negative and 
does not imply the oscillatory behavior displayed in Figs. 1 and 2. 

One can gain information about the individual u and d quark densities by invoking charge 
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FIG. 2: Upper panel: F\. Lower panel: bpib) in transverse position space. The solid curves are 
obtained using [17|] and the dashed curves with 161 ]. 
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FIG. 3: Transverse densities for up u and down d quarks. Each is normalized to unity 

symmetry (invariance under a rotation by n about the z (charge) axis in isospin space) j2l| 
so that the u, d densities in the proton are the same as the d, u densities in the neutron. We 
also neglect the effects of ss [221] or heavier pairs of quarks. In this case p u (b) = p p (b) + 
p n (b)/2, Pd(b) = p p (b) + 2p n (b). The results, shown in Fig. 3, and obtained with either form 
factor parameterization are that the central up quark density is significantly larger than the 
central down quark density. The latter becomes larger at distances away from the center. 
Model independent information about parton distributions has been obtained. In particular, 
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the central density of the neutron is negative. Future experimental measurements of neutron 
electromagnetic form factors could render the present results more precise, or potentially modify 
them considerably. Obtaining a quantitative and intuitive understanding of our results presents 
a challenge to lattice QCD and to builders of phenomenological models. 
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